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Space-time analysis: HBT at SPS and RHIC 
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^ ■ In this review talk, I first summarize the conclusions which can be drawn from a particle 

interferometric analysis of CERN SPS data and then I address related questions which 
are relevant for the upcoming experiments at RHIC. 

> . 

22 ' 1- Introduction 

'NT 

■ The CERN SPS lead beam program has covered for the first time with high statistics all 
O . HBT radius parameters of the full three-dimensional parametrization over a wide region 
Q>^ ', in rapidity and transverse momentum (for data, see Ref. []T|). Theory provides by now 
a detailed data analysis strategy and interpretation (for a review, see Ref. 0) . In this 
talk, I summarize the conclusions which we can draw from the recent CERN SPS data, 
Q ■ and I identify questions which can be settled by the upcoming RHIC experiments. 
2 ■ The main goal of HBT is to determine the geometry and dynamics of the heavy ion 
collision at freeze-out. This information is becoming increasingly important for the dis- 
cussion of other observables. For example: i) In studies of anisotropic flow one 
^ I aims at exploiting the HBT space-time picture to disentangle geometrical and dynamical 
c5 ■ components of anisotropy 13-0]. These cannot be separated on the basis of one-particle 



spectra alone, ii) In the search for event by event multiplicity fluctuations of qualitatively 
new origin (e.g. DCCs), fluctuations due to Bose-Einstein correlations are an important 
background source 0. Their calculation depends strongly on the phase-space density 
attained in the last stage of the collision. The HBT analysis of experimental data now 



provides estimates of this phase space density for a large variety of collision systems []T0 



and work on its consequences for multiplicity fluctuations is in progress [11]. iii) In event 



generator studies, HBT measurements provide constraints on the simulated space-time 



structure |T2|, and the question how microscopic dynamics can account for the observed 
final state emission region becomes of the utmost importance. 

In the light of these rapidly developing applications, one should not forget the funda- 
mental limitations of any particle interferometric analysis. HBT extracts the phase space 
density S{x, K) from a combination of hadronic one-particle and identical two-particle 
spectra of the abundant hadron species: 

e'^ = fd^xS{x,p), (1) 



^ '^^ " Id^xS{x,P,)Jd^yS{y,P,)- 



2 



flow determines ( 



and size of 




Figure 1. Characterization of the pion emission function (^) in terms of Gaussian widths 
R, At], At. The strength of the collective flow field determines the correlation between 
points of highest emissivity and predominant emission direction. The temperature T 
determines the degree of thermal smearing around these points of highest emissivity. 



However: i) S{x, K) characterizes the geometry and dynamics of the collision at freeze-out 
and thus provides only a snapshot of the last stage of the collision. What is characterized is 
the result of a dynamical evolution, not the dynamical evolution itself, ii) HBT measures 
the relative distance between emission points. It does not measure the absolute position 
of emission points. This is of particular importance for the interpretation of the average 
emission time tq, see below, iii) The measurement of C{K,q) gives independent access 
to only three of the four space-time dependencies of S{x,K). This is a consequence of 
detecting on-shell particles. It implies a model-dependent analysis strategy in which a 
model emission function S{x, K) is compared to data. 



2. Final State Geometry and Dynamics of Pb-|-Pb at 158 GeV/c 

The expression for the one-particle and two-particle spectra (|l]) and (H) in terms of the 
emission function S{x, K) form the starting point for various different approaches: i) In 
hydrodynamical simulations of heavy ion collisions ITSlJT^ , they specify the integrations 



over the freeze-out hypersurface. ii) For event generator studies, a discretized formulation 



exists 1 12]. iii) Analytical model studies start from a simple parametrization of the model 
emission function S{x, K) whose model parameters are then extracted from a fit of (|l|) 
and (0) to data. A simple, frequently used model ansatz for the phase space distribution 
S{x, K) in terms of very few, physically intuitive fit parameters, is e.g. 

5(x,P) = tFP-n(x)exp(-^)exp(-^-^-|i^) . (3) 

The physics contained in this ansatz is described elsewhere in much detail @]. As shown 
in the sketch of Fig. |I|, the spatial extension of the pion emission region is characterized 
by Gaussian widths i?, Ai] and At in the transverse, longitudinal and temporal direc- 
tion respectively. The absolute source position in the Minkowski diagram is fixed by the 
time scale Tq, and the pion source emits particles at temperature T. The model allows 
for a directed dynamical component via the flow field u^{x) which is assumed to show 
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Figure 2. Confidence levels for the optimal temperature T and transverse flow r]f val- 
ues from fits of particle emission functions to NA49 one-particle spectra (dashed lines) 
and two-particle spectra (solid lines). The transverse density profile of model emission 
functions was chosen box-shaped (l.h.s.) and Gaussian (r.h.s.). Figure from Ref. [p3 . 



Bjorken scaling in the longitudinal direction superimposed with a transverse flow compo- 
nent sinh?7t(r), ?7t(r) = rjfr/R. To sum up: the phase space density (H) is characterized 
by only 6 physically intuitive fit parameters: T, rjf, R, At], At and tq. 

Studies of many modifications of the parametrization (^ exist. These include e.g. 

in 



the effects of resonance decay contributions [[15 



exist, 
other flow profiles 



temperature 



gradients 17.1^, source opacity and azimuthal asymmetries [P]. Comparisons of the 
particular model (^ with experimental data exist for AGS and SPS [pT] , p2| energies. 
The main argument in favour of a largely model-independent interpretation of the model 
parameters extracted in these studies is that the measured particle spectra (|1]) and (g) 
are mainly sensitive to the various r.m.s. of the emission function and depend only 
weakly on the particular analytical implementation. What matters is mainly the average 
transverse flow and the r.m.s. width of the transverse density profile, not their particular 
r-dependence. A recent analysis [^] of NA49 data indicates the limitations to these 
statements: 
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has compared the model emission function (^ to a model in which 
the Gaussian transverse density distribution of (0) was replaced by a box-shaped one, 
exp (—2^) — ^ ®('" " Rb)- For a comparison of the determined fit parameters of both 
models, one averages the transverse flow over the transverse density profile -ffprofiie('") and 
one recalls that the r.m.s. width for a Gaussian and a box profile differs by a factor 2: 



jj= J rdr tanh(?7t(r))iJprofiie( 



V2 Rg = -^Rb 



Vt{r) = Vf-, 



(4) 
(5) 



The slope of the transverse one-particle spectrum is due to a combination of random 
(thermal) and directed (flow) transverse motion. Different combinations of temperature 
T and transverse flow rjf lead to the same effective blue-shifted temperature. As seen 
in Fig. ^, the parameter combinations favoured by the fit specify a sharp valley in the 
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Figure 3. Comparison of the NA49 Pb+Pb HBT radius parameters at mid rapidity {h^ 
and h~ h~ data are averaged) to results from model emission functions with box-shaped 
(l.h.s.) and Gaussian (r.h.s.) density profile. Different curves indicate the parameter sets 
shown in Fig. 0. Figure taken from . 



T — ?7/-plane. The contour plots in this figure are similar for a Gaussian and box-shaped 
density distribution, but they are not identical. The reason is that for an r-dependent 
flow field, particle emission from different transverse distances r occurs with different 
average transverse momentum. Modifying the transverse density profile -ffprofiio(^) can 
thus affect the prdependence of the transverse one-particle spectrum. For the low-pt 
region [pt < 1-5 GeV) measured by NA49, these model-dependent variations are seen to 
be small. However, the high pt tails of the measured spectrum test in these analytical 
models the large r-tails of the distribution i/profiic('") which is almost unconstrained by 
other observables. An application of these models to higher pt (see e.g. Ref. ||2^) may 
hence be expected to emphasize the model-specific r-dependence of -ffprofiie('") and rit{r). 
The model-dependence increases with pt - this is an important message for similar studies 
at RHIC. 

HBT correlation radii allow to disentangle the T — r]f ambiguity of the fit to the one- 
particle spectrum. This is seen from the contour plot of Fig. |^. The contours for the 90, 50 
and 10 % confidence levels are still very wide due to current experimental uncertainties, 
but they clearly show a different T — rjf correlation. The fit favours a box-shaped density 
profile over a Gaussian one. The parameter sets indicated by bl, b2, b3 and gl, g2, g3 in 
Fig. 1^ are used as input for the fits to the HBT radii in Fig. ^. In the T — rjf plane, the 
fit is driven by the slope of the side and out radius parameters, and the preferred values 
thus depend strongly on taking all contributions to this slope properly into account. All 
existing data analyses p0|^3| favour relatively low temperatures and large transverse 
flows. 

The values of the model parameters used for the fits of Fig. ^ are summarized in the 
table at the end of this section. The box-shaped and Gaussian model emission functions 
differ in several details. Especially, a box-shaped density profile reproduces the data with 
better x^- A similar trend was also pointed out in a recent study of deuteron coalescence 
models [^. But irrespective of these finer tendencies in the data, the main message is 
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unaffected by details of the analytical form of the emission function: we observe a system 
with strong transverse flow which has expanded up to freeze-out to twice its initial size, 
and which has cooled down substantially during this expansion. 

It is an important problem to understand the dynamical origin of these source char- 
acteristics. One may e.g. question the dynamical consistency of the presented data on 
the ground that the average transverse velocity v± and the emission time Tq seem barely 
enough to account for an expansion to twice the initial size. However, HBT does not 
measure absolute positions, and the model parameter tq thus provides only a lower bound 
for the total time after impact 0. Also, very little is known about the question how 
fast the final transverse flow of a hadronic system can be built up during the collision. 
At present, one cannot even rule out the point of view that without novel dynamical 
assumptions, none of the existing event generators can account for the observed strong 
transverse expansion solely with final state interactions |^ . 



A constructive first step to clarify the dynamical origin of HBT fit parameters is the 
comparison of HBT radii from event generators and experiment combined with a compar- 
ison of the simulated space-time structure to the model (0). This may help to understand 
how a dynamical model in which the time after impact is known can account for the 
measured data. 
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3. Phase space density 

Once the emission function 5'(x, K) is determined, the particle phase space density 
/(x, p,t) is obtained by integrating S{x,K) over all particles emitted up to the time 
t, 0]. Of particular interest is the spatially averaged phase space density (/)(K) which is 
directly related to the measured particle spectra p7||28|,|l0|: 



^■''^ ^ Jd^xf{^,K,t>tf) dY M^dM^d<!> Ves{K^,Y)' 

This is a measure for the number of particles in a momentum bin dY dM^ which 
are contained inside the spatial volume Ves_{K^,Y). The spatial volume Vea{K^,Y) is 
determined by the measured HBT radius parameters: 

Kfr(i^x, Y) = ^^^-^^ [Rs ^RIR1-{RI?) . (7) 
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Figure 4. Spatially average phase space density in narrow rapidity windows as a func- 
tion of the transverse momentum. Dashed lines denote Bose-Einstein distributions at 
temperatures T = 80, 120 and 180 MeV. Figure from (l.h.s.) and (r.h.s.). 



The A intercept parameter in 
contributions 
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) ensures that pions from long-lived resonance decay 
which are emitted outside the collision region, do not contribute to the 



average phase space density (/)(K) 



With the help of (j^), D. Miskowiec found in the E877 data a factor 2 decrease of the 
phase space density in the very forward rapidity region. A comprehensive analysis [|10 



of various different collision systems reveals now that such strong rapidity variations of 
(/)(K) occur only very close to the kinematically allowed boundaries. As seen from 
Fig. ^ except for the projectile rapidity region, (/)(K) is almost universal. The variation 
of (/)(K) is small compared to the order of magnitude change in the dN~ /dY between 
the different collision systems (22 for S-S vs. 185 for Pb-Pb). The transverse momentum 
dependence of (/)(K) is in rough agreement with a Bose-Einstein distribution. Deviations 
from this show up in a weaker transverse momentum dependence of (/)(K) and can be 
accounted for by transverse flow. This is seen on the r.h.s. of Fig. ^ where the average 
phase space density extracted from NA49 mid rapidity data is compared to results from 
the optimal model parameter sets 'bl' and 'gl' indicated in Fig. 0. 



4. Extrapolation of the AGS and SPS systematics to RHIC 



All observations indicate that the decoupling of pions occurs at an approximately con- 
stant phase space density irrespective of the collision system. This implies that the freeze- 
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out volume scales with particle multiplicity 



dYMi dM, rf$ 



(8) 



The qualitatively new aspect of Fig. ^ is to establish equation @ differentially, i.e. for 
narrow bins in transverse momentum and rapidity, i^^-integrated versions of are 
well-known. NA35 observed e.g. for different collision systems a linear increase of the 
averaged collision volume with dN~ /dY up to dN~ /dV 



60 pg. TheKj 



averaged 
180, and 



NA49 data for central rapidity follows the same systematics up to dN^ /dY 
RHIC will allow to test this systematics even further. 

Eq. (§) can be used for extrapolation into this new high multiplicity regime accessible 
to RHIC: The most naive RHIC-prediction is obtained by turning Ves oc dN/dY into a 
simple multiplicity dependence for the side radius parameter which is typically the least 
affected by the dynamics, Rs oc {dN/dY^^^. This leads to i?, ^ 9 fm for ^ = 600 
and i?s ^ 11 fm for = 1200. Depending on the dynamics, however, even if (||) holds, 
the transverse radius parameters may be larger or smaller. One such example, based 
on a hydrodynamical picture, assumes that the transverse flow gradient does not change 
between heavy ion collisions at SPS and RHIC. From this picture, one expects changes 
in the emission time Tq, transverse radius R, and transverse flow rjf 
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The resulting HBT radius parameters indicate a collision system which is substantially 
more extended along the beam axis at the expense of a relatively small increase in the 
transverse extensions. 
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(10) 



Most scenarios favour a sudden bulk freeze-out of all pions which is reflected in small 
values for R^ — R^. However, hydrodynamical scenarios with a particularly soft equation 
of state lead to large differences between the two transverse parameters R^ R^. 
In these scenarios, only ^ 50 percent of the difference R^ — R^ stems from an extended 
lifetime of the system. An equally important contribution to R^ — R^ comes from a strong 
negative x-t correlation \T^\. 



5. Multiparticle Bose-Einstein symmetrization 

Standard calculations of the two-particle correlator C {K, q) are based on approximating 
A^-particle symmetrized flnal state wavefunction by a sum of two-particle symmetrized 
contributions. For small phase space densities, this is a valid approximation. However, 
if the phase space density is sufficiently high, multiparticle symmetrization effects cannot 
be ignored. Recently, progress was made in quantifying the decisive words "sufficiently 
high" of this last sentence in simple model studies: 

The well-known combinatorial problem is that one obtains from an A^-particle sym- 
metrized wavefunction A^! contributions to the measurable particle spectra. Even for mod- 
erate event multiplicities, brute force numerical calculation do not work (100! ~ 9*10^^^). 
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However, for pion emitting sources defined by the emission function S{x, K) and the event 
multiphcity A^, Pratt's algorithm M provides a substantial simplification: 



Input: S(x,K) 
multiplicity N 



Pratt terms: 
G„(Pi.Pj) 



ij,n 



w(n), 
: [1,N] 



Output: Spectra = 
Z (n (Pratt terms)) 

partitions 



Here, the building blocks G„(pj), C„ and w{n) are defined in terms of n-dimensional 
integrals over S{x, K), n < N. All measured spectra reduce to sums over products of these 
building blocks. While the technicalities of this approach cannot be explained shortly (see 
Ref. [0 for a review), the important message is that the sums involved in the final answer 
run over the number of all partitions of different contributions, rather than over A^!. 
Note that Partitions[N=100] =2 * 10^, and this makes a numerical calculation feasible. 

The remaining computational task is the calculation of the n-dimensional integrals Gn, 
Cn and w{n). For Gaussian emission functions, this can be done analytically [p2i^ . 
Fig. ^ shows the results of such a model study and confirms the observations of Zajc . 
Pions like to sit close together in momentum space and this leads to a steepening of the 
one-particle spectrum. Also, they like to sit close together in configuration space and 
hence the two-particle correlator gets wider. As can be seen from Fig. |^, these effects are 
small but visible for realistic phase space densities below 1 h^. 

However, existing calculations like Fig. |^ cannot be compared to experiment. The 



problem is that despite some recent progress [35|, there is no realistic description of muL 



tiparticle final state Coulomb interactions. The HBT effect and Coulomb repulsion occur 
both on the same relative momentum scale and counteract each other. One may hence 
expect that their multiparticle contributions cancel to a large extent. Thus, present calcu- 
lations can be expected to provide a loose upper bound on multiparticle symmetrization 
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Figure 5. Modifications of one-particle spectrum (l.h.s.) and two-particle correlator 
(r.h.s.) due to multiparticle symmetrization effects. With increasing phase space den- 
sity, the one-particle spectrum gets steeper and the two-particle spectrum gets wider. 
Figure taken from p3| . 
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effects. The important message for RHIC is that if the so far universal value of (/)(K) is 
confirmed at RHIC, multiparticle effects to the momentum spectra are of minor impor- 
tance. 



6. Multiplicity fluctuations 

There is one important observable where the calculation of multiparticle symmetrization 
effects allows for a comparison with experimental data and that is the modification of the 
charged to neutral multiplicity distribution due to Bose-Einstein statistics. A typical 
Bose-Einstein modified distribution, proposed originally for the description of anomalous 
cosmic ray (CENTAURO) events, is 



1 iV! w{no) w{n+) w{n^] 

P{no,n+,n-) = 



Z 3^ uqI n+! 71- 



(11) 



Here, the total number of neutral and charged pions is kept constant, but one assumes 
that the multiplicity distribution is arranged according to BI3-statistics via the channel 
7r°7r° — vr+vr". Z is a normalization, and w{n) denotes the multiparticle symmetrization 
weights which include higher order correlations. According to Pratt's formalism, they are 
given by a sum over all partitions (m, lm)n of the building blocks Cm 

w{n) = — ^ — (12) 

(m,/m)n Hm "^"^ (<m-) 

These expressions simplify substantially for Gaussian models of the emission function 
S{x, K). One observes that the leading contribution to the building blocks Cm is a simple 



polynomical in the phase space volume e |33 



h 



3 



Cm^e""-', where e = -—^ , (13) 



and this reduces the complicated combinatorics of (^) to a simple polynomial 

n 

w{n) = + t{k - I)) . (14) 



k=l 



Based on equations (|TTp - ([T^), P. Steinberg |TI| has studied the charged to neutral parti- 
cle ratio measured by WA98. He starts from the observation that the distribution of this 
ratio is approximately 15 % broader than that obtained in a VENUS 4.12 event generator 
simulation. Clearly, this discrepancy can have many reasons. The strong working hy- 



pothesis in Ref. |11] is that the event generator describes the physics correctly but misses 
multiparticle symmetrization effects (as all Monte Carlo models do). As a remedy of the 
latter, the VENUS 4.12 multiplicity distribution was modified via (|ll|). P. Steinberg finds 
that for phase space densities consistent with those shown in Fig. ^, multiparticle corre- 
lations can account for a substantial part of the observed discrepancy (approx. 8-10 %). 
Irrespective of the real origin of the discrepancy, this is a remarkable result: even for the 
relatively small measured phase-space densities, Bose-Einstein statistics leads to visible 
changes in multiplicity distributions. This also indicates that Bose-Einstein correlations 
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are an important background source in searches of qualitatively new physics (e.g. DCCs) 
in multiplicity distributions. 

I thank U. Heinz, J.G. Cramer, D. Ferenc and B. Tomasik with whom I had the privilege 
to collaborate. This work was sponsored by DOE Contract No. De-FG-02-92ER-40764. 
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